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Abstract. A simple analytically solvable model exhibiting a 
1// spectrum in an arbitrarily wide frequency range was recently 
proposed by Kaulakys and Meskauskas (KM). Signals consisting 
of a sequence of pulses show that inherent origin of the 1// noise 
is Brownian fluctuations of the average intervent time between 
subsequent pulses of the pulse sequence. We generalize the KM 
model to reproduce the variety of self-affine time series exhibiting 
power spectral density S{f) scaled as power of their frequency 
/. Numerical calculations with the generalized discrete model 
(GDM) reproduce power spectral density S{f) scaled as power 
of frequency 1//^ for various values of /?, including /5 = 1/2 for 
applications in financial markets. The particular applications of 
the model proposed are related with financial time series of share 
volume traded. 



1 Introduction 



Physicists have recently begun doing research in finance with wide apphcation 
of models earlier introduce in physics ^ H . The statistical properties of 
financial time series are attracting experts of statistical physics, fiuctuations 
analysis, dynamical chaos and others. Papers on finance are appearing with 
some frequency in physics journals. A new movement called econophysics 
has been established [0. On the other hand, several empirical studies have 
determined scale-invariant behavior of both the long range correlations of 
price volatility and share volume traded in financial markets p, ^ |^, ^. 
Mandelbrot introduced the concept of fractals in terms of statistical self 
similarity [0 and using the context of self-affine time series extended the 
concept to time series The basic definition of self-affine time series is 
that the power spectral density of the time series has a power-law dependence 
on frequency S{f) = f~^. Universality of 1// noise, when /? = 1, has led 
to speculations that there might exist some generic mechanism underlying 
production of so general statistical properties. Such concepts as fractional 
Brownian motion provide a formal procedure how to produce self-affine time 
series, but can't serve as generic mechanism. We will generalize the simple 
model introduced by Kaulakys and Meskauskas (KM) to generate time 
series in the range < /9 < 2 with particular interest in financial time series. 

Long range correlations in time series I{t) are quantified by autocovari- 
ance (autocorrelation) function C(s): 

c{s) = ci-s) = {fiy mnt + s)dt'j , (1) 

with Wiwner-Khintchine relation to power spectral density S{f) defined as: 
S{f)=Yiml% ! I{t)e-'^^f'dt )=4 1im I C{s) cos{27ifs)ds, (2) 

T^oo \ T Jo / Jo 

where T denotes the considered time interval of time series I{t) . The KM 
model is based on the time series generated as: 

Ht) = Yl "Ik^it - tk) (3) 

k 

where is a contribution to the signal of one pulse and noise is due to the 
correlations between the occurrence times of 5 type pulses. This model 
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corresponds to the flow of point objects: pliotons, electrons, cars, trades in 
financial markets and so on. The simplest version of the model consist of 
sequence of transit times tk described by recurrence equations: 

tk = tk-l + Tk, 

Tk = n~i - l{rk-i + aSk. (4) 

Here the recurrence time = tk—t^-i fluctuates due to the external random 
perturbation of the system by sequence of uncorrelated normally distributed 
random variables {sk} with zero expectation and unit variance, where a 
denotes the standard deviation of the white noise, 7 <^ 1 is the recurrence 
time r relacsation rate to the some average value r. Note that r follows an 
autoregressive process AR(1). 

This model containing only one relaaxation time 7"^ can for sufficiently 
small parameters 7 produce an exact 1// - like spectrum in wide range of 
frequence 0: 

SU) = ^, /i</<min(/2,^), (5) 

'T f 

where an is a dimensionless constant — the Hooge parameter: 

2^e-^', K=^. (6) 

and /i = ^^/yira, h = 2^^l^/7:a, U = (2vrr)-i. 

Here we present generalized KM model with particular interest to repro- 
duce the statistical properties of share volume traded in financial markets. 
We will generate discrete share volume time series with recurrence time of 
distinct trades described as AR(2) process. This particular application of the 
model proposed by B. Kaulakys and T. Meskauskas will exhibit univesality 
of the mechanism underlying production of an noise. 

2 Model definition 

In the numerical data analysis we usually deal with discrete sets of data. 
There is the first need to modify the model to the discrete one. Let us 
introduce a new conventional time scale defined with a time interval r^. 
Integrating continuous signal I{t) in the subsequent intervals of length 
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we will get discrete time series (let us call it volume Vr) 



Vr 



I I{t)dt = ^Qk, tr = TTd. (7) 



Consequently, the discrete power spectral density can be expressed as: 



S{fs 




r=l 



(8) 



where fs = = Tan. 

New discrete series are equivalent to the initial sequence of S functions, 
when Td whit the same values of 7 and a. Numerical results of S{fs) 
calculated with various parameters are shown in Fig. |^. Multiple numerical 
calculations confirm full correspondence of the discrete model (DM) defined 
here with earlier introduced by B. Kaulakys and T. Meskauskas 0. 

The model can produce the 1//- like spectrum in an arbitrarily wide 
range of frequencies /i < / < /2 , /t and is free from unphysical divergence 
of the spectrum at / — 0; for f < fi- 

^(^) = " l + aV4rV 

Due to the long memory random process, defining transit time sequence 
tk, the model describes long time correlations quantified in power spectral 
density S{f) ^1//. This model may be also generalized for non-Gaussian 
distribution of the periods r^. Then 

S{f) = 2r-V(0)//, (10) 

where iplr) is the distribution density of periods r^. This makes the model 
applicable to the wide variety of stochastic processes, which have well defined 
distribution function ^(r) in the vicinity of = 0. Numerical calculations 
confirm that, when ip{0) = 0, the dependance of the power spectral density 
on frequency appears as S{f) ^ l/P^"^ |T0| . 

One more generalization of the model is needed for applications in finan- 
cial markets and other self-affine time series with a power-law dependence 
on frequency S{f) = f^^ with /? ~ 1/2. We will strengthen high frequencies 
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S{f) S(f) 




Figure 1: Power spectral density versus frequency calculated from the model 
described by Eqs. (3), (4), (7), (8). The main parameters defining the 
model are: r, a, 7. Time scale determines the highest frequency f < ^ 

under consideration. The sinuous curves represent the results of numerical 
simulations averaged over five realizations, and the straight lines represent 
the analytical spectra described by Eqs. (5) and (6). 
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Figure 2: Volume Vr versus the number r of time interval and power spec- 
tral density S{f) versus frequency /, calculated from the generalized discrete 
model (GDM) with a = 0.0002, r = 0.01, a = 0.95,7 = 0.0002, = 1. S 
is averaged over five realizations. The straight lines represent the fractional 
power law S{f) = 2* 10^/ f^^^. 



and will account for positive playback of Tk increment by adding the term 
aArfc_i = a(rfc_i — to the recurrent expression: 

Tk = Tk-i + a(rfe_i - Tk-2) - lirk-i - r ) + aSk- (11) 

Note that this new term changes autoregressive process AR(1) to the higher 
one AR(2). Multiple numerical calculations with the generalized discrete 
model (GDM) exhibit dependance of /? on a and other parameters of the 
model: cr, r, 7. We demonstrate an example of numerical calculation with 
GDM in Fig. |^, which exhibits clear fractional power law with (3 = 1/2 oi 
the power spectral density. 



3 Application to the financial market 

An important quantity that characterizes the dynamics of price movement 
in financial markets is the number of shares Vr (share volume) traded in 
a time interval rT^i < t < {r + l)T(f. The statement "It takes volume to 
move stock prices" accumulates very general idea that statistical properties 
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of financial markets are enclosed in time series of share volume. Very direct 
confirmation of this statement and quantitative investigation of the largest 
1000 stocks in three major US stock markets was recently presented in 0. 
This work provides an evidence that long range correlations in share volume 
and price volatility are largely due to those of the number of trades Nr in 
time interval r^. Close correlation between A^^ and Vr is imposed by the 
relation: 



and very weak correlation in time sequence of share volume per transaction g^. 
These results suggest us to apply GDM as a model for the time series of share 
volume traded in financial markets , with a simple assumption that in the 
first approach the average (qi) can be included into the normalization factor. 
This mean that simple relation Vr = Nr {qt) enables us to make comparison 
of GDM to the variety of real market data. In Fig. ^. we demonstrate 
comparison of Lithuanian Stock Market data with numerical results from 
GDM. The volume of shares included in the index LITIN is normalized to 
average of 100 trades per = Iday . Note that despite the model simplicity it 
serves as market data generator and reproduces the main statistical property 
of the system, i.e., the power law dependence (/3 = 1/2) on the frequency of 
the power spectral density. 

4 Conclusion 

The empirical evidence provided by Gopikrishnan a/ Q, that the num- 
ber of transactions A^^ in the subsequent time intervals define long range 
correlations of share volume traded, enabled us to apply simple model of a 
1/f noise and to reproduce long-range correlations of the share volume 
traded in financial markets. We generalized the KM model by integrating 
the sequence of pulses in a conventional time interval and replacing the 
recurrent time Tk stochastic process from AR(1) to the AR(2) one. Numeri- 
cal calculations with the generalized discrete model (GDM) reproduce power 
spectral density S{f) scaled as power of frequency 1/ for various values of 
f3, including /3 = 1/2 for applications in financial markets. Further investi- 
gation of the model with its possible applications in financial time series is 
in progress. 




(12) 



i=l 
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Figure 3: (LSE) The volume of shares , included in the index LITIN 
of Lithuanian Stock Exchange, normalized to an average of 100 trades per 

Tfi = Iday versus the number of day traded and corresponding power spectral 
density S{f) versus frequency / calculated from Fast Fourier Transform of 
discrete data. (GDM) Volume Vr versus the number r of the time interval 
Td — 1 and power spectral density «S'(/) versus frequency /, calculated from 
the generalized discrete model (GDM) with a — 0.00015, r = 0.01, a — 
0.75, 7 = 0.0007, Td — 1. The straight lines approximating the power spectral 
density curves represent the fractional power law S{f) — 700//^^^. 
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